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Martensites are metastable phases, possessing a characteristic morphology, usually formed during 
a fast quench accross a structural transition. We attempt to understand these morphological features 
using a coarsegrained free energy functional T[e; $] which contains, in addition to the usual strain 
fields eij (the " order parameter" for the transition), the "vacancy" field <f) which arises due to 
the geometric mismatch at a parent-product interface. The relaxation of this mismatch is slow 
compared to typical front propagation times and hence <f) is essentially frozen in the reference frame 
of the growing martensite front. Minimisation of T then automatically yeilds typical martensite 
morphologies. We demonstrate this in two dimensions for the square to rhombus transformation 
and obtain internally twinned martensites, which grow as thin strips, for "hard" martensites (e.g., 
Fe-based alloys) or with a 'single-interface', for "soft" martensites (e.g., In-Tl alloys). 

a" 

I. INTRODUCTION AND PHENOMENOLOGY 

> : 

Martensites, a generic term applied to diffusionlcss structural transformations producing long-lived metastable phases 
with characteristic morphologies, are observed in a wide variety of systems 0^], ranging from metals, alloys and 
ceramics to organic crystals and crystalline membranes || . Although a substantial amount of empirical knowledge 
has been acquired over the years, a systematic theoretical study of the morphology and kinetics of martensites |^,^) is 
\^ " lacking. Analysis is made difficult by the fact that the formation of martensite is necessarily a non-equilibrium process. 
ON ■ A theory of the martensitic transformation should, therefore, be able to describe the system for all times, t, and strictly 
equlibrium (t — > oo) descriptions are invalid. In this paper, we derive from microscopies a free energy functional 
capable of describing a droplet of growing martensite at intermediate times. We obtain several, empirically observed, 
p] morphological features of Martensites by minimising the free energy functional subject to appropiate constraints. In 
the rest of this section we review relevant aspects of Martensite phenomenology. We derive the droplet free energy 
functional in Section 2. In Section 3 we use it to investigate the morphology of two-dimensional martensites. Lastly, 
in Section 4 we discuss our results and conclusions. 

As a typical example of a Martensitic transformation consider the case of steel (an alloy of iron, Fe and carbon, C). 
Pure Fe at atmospheric pressures undergoes an equilibrium structural phase transition from the austenite (7) phase, 
which has an FCC structure to a BCC, ferrite (a) phase, at a temperature T c = 910°C. A fast quench from T > T c to 
T < T rns (martensite— start temperature) < T c , produces, however, a rapidly transformed metastable phase called the 
martensite. On nucleation, martensite "plates" grow with a constant front velocity (v ^10 5 cms~ 1 in Fe) || which 
is comparable to typical acoustic speeds in solids. The plates consist of an array of twinned or slipped BCC crystals, 
having well defined orientational relations with each other and with the parent matrix. The Miller indices of the habit 
plane (the equatorial plane of the plates), are fixed relative to the parent. The morphology of the plates depends on 
the elastic moduli of the solid, ranging from "lens-shaped" discs (typically ~ l^im) with a high aspect ratio in the 
case of hard Fe-based alloys to a region bounded by a single interface in the case of the softer In-Tl alloys. If the 
clastic energy barriers are large, the plates stop on colliding with other plates. The resulting pattern displays a scale 
invariant size distribution of the plates j7j . 

Several empirical facts suggest that martensites are always metastable with respect to the ferrite. Firstly, the 
detailed morphology of the plates is a sensitive function of the quench history. Secondly, mechanical properties of 
martensites are dependent on aging and tempering — metallurgical processes which consist of holding steel at elevated 
temperatures for specific periods of time to obtain required properties. Lastly, martensites are never formed when 
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steel is cooled sufficiently slowly, thus there is an intimate relation between the kinetics of cooling and the morphology 
— slow cooling nucleates a spherical droplet of the ferrite, while rapid quenching produces a martensite with a well- 
defined morphology. Any theory of the Martensitic transformation should therefore include variables that change 
slowly over time scales corresponding to the nucleation front propagation. Such typical "slow modes" of a solid || 
undergoing a first-order structural transformation are the displacement field u(r,i) (goldstone mode), the mass and 
momentum densities p(r, t) and g(r, t) (conservation laws) and the strain tensor Cij(r, t) (the symmetry breaking order 
parameter). Imagine, however, a grain of the product nucleating within a parent matrix at time t = 0. 
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Fig. 1 Inclusion of a product (rhombus) in a parent matrix, showing the discontinuity in u across 
the parent-product interface at t = 0. 

Fig. 2 Variational ansatz for a polycrystal of length L and width W showing N — 4 twins along 
the x axis. The regions marked + and — are related to each other by a mirror reflection at the twin 
interface. The interfacial widths rj, £ and £ are also shown. 

It is clear (Fig. 1), that an atomic mismatch is generated at the parent-product interface || as soon as the nucleus 
is formed. This mismatch a ppe ars as a discontinuity in the normal component of the displacement field across 
the parent-product interface |l2f| (a consequence of the fact that the 7 — > a structural transformation involves a 
spontaneous breaking of a discrete symmetry) , Au • n = <fi, and leads to a local compression or dilation of the atomic 
environment at the interface. Since in martensites, the transformation front velocity, v, is large compared to atomic 
diffusion times, </>(r, i) appears as a new slow mode (vacancy field). Once formed, this new field, relaxes diffusively 
over a time scale >> L/v, where L is the typical size of the grains. A Langevin description of the dynamics of 
the front, requires a free-energy functional which describes all intermediate configurations between an austenite and 
a ferrite. Thus the free-energy functional should include configurations corresponding to crystalline inclusions of the 
a phase in the 7 phase, where the vacancy field <f> ap pea rs explicitly. The usual elastic free-energy functional of a 



solid, F e i, has to be augmented by an interfacial term |12|, Fi nt , describing the parent-product interface. It must be 



stressed that these droplet configurations do not affect equilibrium behaviour. Thus the equilibrium ferrite will be 
obtained by minimising F e \. The dynamics, however, is crucially altered by the inclusion of such fluctuations. In the 
next section we derive, starting from a microscopic density functional theory, a free energy functional which includes 
such droplet fluctuations. 

II. THE DROPLET FREE ENERGY FUNCTIONAL 

Consider a droplet of the product phase of size L within a parent matrix bounded by an interface of thickness £. The 
presence of the droplet results in a strained parent crystal such that the strain varies from zero far away from the 
droplet to a limiting value at its center. The free-energy functional JlO) ] of this strained parent is given by, 

c(|R-R' + u(R)-u(R')|), (1) 
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where {R} represent the lattice vectors of the parent crystal, u(R) are the displacement fields (deviation of each atom 
from the (perfect) lattice positions of the parent), and c(|r — r'Q is the direct correlation function [jllj of the liquid at 
freezing. The above expression is exact at T = ; corrections are of the order of the r.m.s. fluctuations of the atoms 
about their perfect lattice positions, which are small in the solid phase. 

Since the range of c is of order £, the sum over {R, R'} in Eq. |l|, can be split into three parts — (i) RAR' € parent , 
(ii) RAR' € product and (iii) R € parent and R' G product. Expanding u(R') about R in the first two parts, 
leads to the usual bulk elastic free-energy functional F e i. However, as argued above, the assumption of continuity 
breaks down in (iii), and so such an expansion cannot be carried out. In the limit £/L <C 1, the interface can be 
parametrized by a sharp surface T (r) = 0, with R and R' lying infinitesimally close to T = 0. Across this interface, 
the normal component of u is discontinuous ; Hto|r_r/|^ ( u (R-) — U (R-')) 1 n = (u+ — u_) • n = 4>5(T), where n is the 
unit normal to the interface. The discontinuity </>(r) = (p(r)\r=o — Pav)/Pav is the local vacancy field, while p is the 
local density and p av is the density averaged over a unit cell of the parent. In the continuum limit, the contribution 
to T from (iii), to leading order in the discontinuity, reduces to 

^ J dr[(u + -u_).n] 2 ,5(r). (2) 

The coefficient 7 = j dr S(C,)^ 2 c" (r) is clearly the inverse surface compressibility modulus of vacancies, 4>. This is 
expected to be of the same order of magnitude as the other elastic moduli, and thus ranges from 0(1O~ 3 ) for soft 
materials (e.g., In-Tl) to 0(1) for hard materials (e.g., Fe-Ni). Since the strains = (9,tij + djiii)/2 are continuous 
across the interface, we replace the delta function in Eq. ^ by a regulator [1 — exp(— £<9„ejj) 2 ], so that the final 
regulated free-energy functional which incorporates all the slow modes in the problem is, 

T = F el + F mt = F e i + ^ J dr 4? (d n eij f . (3) 

At the intial time (t = 0), the transformed region (product) is simply obtained as a geometrical deformation of the 
parent, Fig. 1, which fixes the initial value of 4>. Having created this discontinuity <f> at the parent-product interface, it 
will diffuse over a time r^. Note the two time scales relevant to our kinetics — the quench rate r _1 = (dT/dt)/T and 
the vacancy relaxation time r^. Accordingly, two extreme dynamical limits, suggest themselves — a "slow quench", 
t > t^, where the <j) field relaxes instantaneously; and a "fast quench", t < t^, where the <fi field is "frozen" in 
the frame of reference of the moving nucleation front and (j) relaxes only once the growth of the front stops [fl3|| . 
Consequences of an infinitely slow quench is well known from standard classical nucleation theory |l4| since the 
<p field relaxes to zero instantaneously and hence does not contribute. In the next section we study infinitely fast 
quenches which produce martensitic morphologies in the simple setting of a two-dimensional first order structural 
transformation. 

III. MORPHOLOGY OF 2D MARTENSITES 



We study the first-order structural transformation from a square (austenite) to a rhombus (ferrite) unit cell in two 
dimensions [ fL5|| . Our results, which can be easily extended to the tetragonal to orthorhombic transition (essentially a 
2-dim square to rectangular [Q), are of relevance to structural transformations in alloys like In-Pb, In-Tl, Mn-Fe. This 
transformation involves a shear+volume deformation, and so the strain order parameter eij has only one nontrivial 
component e3 = (u xy + u yx )/2. 

The bulk elastic free energy has three minima — one corresponding to the undeformed square cell (e3 = 0) and the 
other two corresponding to the two variants of the rhombic cell (e3 ~ ±eo). This immediately leads to the free-energy 
functional F = F e i + Fi n t (in dimensionless variables) , 

F el +F int = J dxdy [ae 3 2 - e 3 A + e 3 6 + {(d x e 3 ) 2 + (d y e 3 ) 2 } + 1 {4>d n e 3 ) 2 ] . (4) 

The three minima of the homogeneous part of F at e3 = (austenite) and e3 = ±eo = ±[(l + \/l ~ 3a)/3] 1 / 2 (ferrite), 
are obtained in the parameter range < a < 1/3. The parameter a is the degree of undercooling and plays the role 
of temperature. 

An infinitely fast quench implies that the <fi field is "frozen" in the frame of reference travelling with the front, 
r — > r — vt and 4>t — > v ■ V</>, where <pt denotes the derivitive with respect to time. Quenching the solid across the 
structural transition, nucleates an inclusion which initially grows as a ferrite. Further growth as a single-domain ferrite 
is discouraged by Fi nt . The growing nucleus, gets around this problem by creating a twin. This costs interfacial energy 
which is however small. We use a variational strain profile to include specifically such twinned domains. For a start, 
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let us consider a rectangular nucleus of length L (along the x-axis) and width W (along the y-axis). The nucleus is 
divided into N twins (Fig. 2) with alternating between eo and — eg. The N— 1 twin interfaces, all of thickness rj, can 
be parametrized by (x), which for the i-th. interface takes values — eo (ferrite) at (i— l)L/N + rj/2 < x < iL/N — rj/2 
and eo (twin ferrite) at iL/N + rj/2 < x < (i + 1)L/N — rj/2, with a linear interpolation in between. The thickness 
of the austcnitc-ferritc interface along the x-axis (y-axis) is £ (£) within which (x) interpolates linearly between 
and eo (or — eo). The free-energy of a nucleus with N twins is 



F(L, W) = [ AFLW + {21 - AF)(L^ - 
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The interfacial widths are taken as variational parameters. The minimised widths £ and 77 are independent of the 
number of domains N, while the width £ scales as L/N. The free energy, once minimized with respect to the widths, 
is a function of L, W and N. As is aparent from Eq.[|, T goes through a minimum as N is varied, for every value 
of L and W. This is a consequence of the competition between the interfacial energy at the twin interface (which 
disfavours twinning) and the interfacial energy at the parent -product interface which favours twinning since it reduces 
<f>. Minimising the free energy with respect to the number of twins TV immediately gives us, L/N — W a (Fig 3). The 
exponent a ~ \ for small W, for large W there are deviations. The exponent a is empirically known Q to lie between 
.4 and .5. A relation similar to ours has been obtained (with a = ^) from equilibrium continuum elastic theory [Q. 
It should be noted here, however, that the statement we make is stronger — we imply that L/N = W a for all times 
during the growth of a martensite plate as a consequence of demanding local equilibrium. Our prediction for L/N 
can be verified from in situ transmission electron microscopic studies of growing martensite fronts. 

The free energy minimised with respect to all the "fast" variables (viz. the interfacial widths and N) is plotted 
in Fig. 4. One observes that for 7 = 1.0 (Fig 4), the free energy, for all choices of L, shows a minimum in W for 
W < L, so that the system will always prefer to have thin rectangular strips reminiscent of acicular martensites seen 
in Fe-Ni or Fe-C systems. For 7 = .001 there is no such minimum and the martensitic region would like to grow in 
both directions. 
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Fig. 3 L/N plotted against W 1/2 for L = 20000 and 7 = .001. 

Fig. 4 Contour plots of the droplet free energy minimized with respect to the interfacial widths and 
iV in the L — W plane for 7 = 1.0. Note that there is a minimum in W for all choices of L (for large 
L and W) in this case. 



Once growth in one of the directions traverses the entire extent of the sample, growth will proceed thereafter only 
along a single interface. Such single interface growth is seen in "soft" In-Tl systems. The critical nucleus is obtained 
by locating the saddle point in each of these free energy surfaces. Using the quantity k = (££/ 'LW) 1 ! 2 as a measure 
of the sharpness of the interface (indeed one can talk of a definite droplet shape only if the size of the droplet is larger 
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than this measure), we find that the critical nucleus is diffuse for both 7 = 1.0 (k = .4) and 7 = .001 (k = 1.43). 
However, as the size of the droplet increases n rapidly reduces to zero in both cases. 

IV. SUMMARY AND DISCUSSION 

In summary, we have derived from microscopies a free energy functional ^[cij, <j>] capable of describing droplets of a 
martensitic product enclosed within a parent austenitic matrix. Essential to this description is a new field (f>(r, t) viz. 
the geometric mismatch between the parent and the product at the austenite-martensite interface. This mismatch 
relaxes diffusively. For martensites, the typical front velocities are empirically known to be very high so that <f>(r, t) 
fails to relax within the timescale of front motion. We have therefore minimised our droplet free energy functional, 
using a variational strain profile, subject to the constraint <j)(r,t) = <p(r, t = 0). We show that typical martensite 
morphologies drop out naturally from such a procedure. The condition of local equilibrium yields experimentally 
verifiable quantitative morphological features. In this context, it is worth noting that in practice, real martensites in 
"hard" systems (eg. Fe-Ni) do not occur as rectangular strips but as lens shaped objects. The following argument will 
show that a droplet in the form of "convex-lens" has, in fact, lower energy. The interfacial energy at either end of the 
rectangular strip (where the interface is along the y direction) is 0( ( W(L) ) = 0(W 2 ), where W(L) is the width of 
the strip at x = L. This can be reduced further by making W(L) (and H^(0), the width of the strip at x = 0) as small 
as possible. The "lens" shaped structure thus obtained is consistent with structures seen in acicular martensites, like 
Fe-Ni and Fe-C. 
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